Autophagy Markers Are Altered in Alzheimer's Disease, Dementia with Lewy Bodies and Frontotemporal Dementia.

The accumulation of protein aggregates defines distinct, yet overlapping pathologies such as Alzheimer's disease (AD), dementia with Lewy bodies (DLB), and frontotemporal dementia (FTD). In this study, we investigated ATG5, UBQLN2, ULK1, and LC3 concentrations in 66 brain specimens and 120 plasma samples from AD, DLB, FTD, and control subjects (CTRL). Protein concentration was measured with ELISA kits in temporal, frontal, and occipital cortex specimens of 32 AD, 10 DLB, 10 FTD, and 14 CTRL, and in plasma samples of 30 AD, 30 DLB, 30 FTD, and 30 CTRL. We found alterations in ATG5, UBQLN2, ULK1, and LC3 levels in patients; ATG5 and UBQLN2 levels were decreased in both brain specimens and plasma samples of patients compared to those of the CTRL, while LC3 levels were increased in the frontal cortex of DLB and FTD patients. In this study, we demonstrate alterations in different steps related to ATG5, UBQLN2, and LC3 autophagy pathways in DLB and FTD patients. Molecular alterations in the autophagic processes could play a role in a shared pathway involved in the pathogenesis of neurodegeneration, supporting the hypothesis of a common molecular mechanism underlying major neurodegenerative dementias and suggesting different potential therapeutic targets in the autophagy pathway for these disorders.

Archival skin biopsy specimens as a tool for miRNA-based diagnosis: Technical and post-analytical considerations.

Archived specimens, taken by standardized procedures in clinical practice, represent a valuable resource in translational medicine. Their use in retrospective molecular-based studies could provide disease and therapy predictors. Microfluidic array is a user-friendly and cost-effective method allowing profiling of hundreds of microRNAs (miRNAs) from a low amount of RNA. However, even though tissue miRNAs may include potentially robust biomarkers, non-uniformed post-analytical pipelines could hinder translation into clinics. In this study, epidermal RNA from archival skin biopsy specimens was isolated from patients with peripheral neuropathy and healthy individuals. Unbiased miRNA profiling was performed using RT-qPCR-based microfluidic array. We demonstrated that RNA obtained from archival tissue is appropriate for miRNA profiling, providing evidence that different practices in threshold selection could significantly influence the final results. We showed the utility of software-based quality control for amplification curves. We revealed that selection of the most stable reference and the calculation of geometric mean are suitable when utilizing microfluidic arrays without known references. By applying appropriate post-analytical settings, we obtained miRNA profile of human epidermis associated with biological processes and a list of suitable references. Our results, which outline technical and post-analytical considerations, support the broad use of archived specimens for miRNA analysis to unravel disease-specific molecular signatures.

Genetic Profiling of Sodium Channels in Diabetic Painful and Painless and Idiopathic Painful and Painless Neuropathies.

Neuropathic pain is a frequent feature of diabetic peripheral neuropathy (DPN) and small fiber neuropathy (SFN). Resolving the genetic architecture of these painful neuropathies will lead to better disease management strategies, counselling and intervention. Our aims were to profile ten sodium channel genes (SCG) expressed in a nociceptive pathway in painful and painless DPN and painful and painless SFN patients, and to provide a perspective for clinicians who assess patients with painful peripheral neuropathy. Between June 2014 and September 2016, 1125 patients with painful-DPN (n = 237), painless-DPN (n = 309), painful-SFN (n = 547) and painless-SFN (n = 32), recruited in four different centers, were analyzed for SCN3A, SCN7A-SCN11A and SCN1B-SCN4B variants by single molecule Molecular inversion probes-Next Generation Sequence. Patients were grouped based on phenotype and the presence of SCG variants. Screening of SCN3A, SCN7A-SCN11A, and SCN1B-SCN4B revealed 125 different (potential) pathogenic variants in 194 patients (17.2%, n = 194/1125). A potential pathogenic variant was present in 18.1% (n = 142/784) of painful neuropathy patients vs. 15.2% (n = 52/341) of painless neuropathy patients (17.3% (n = 41/237) for painful-DPN patients, 14.9% (n = 46/309) for painless-DPN patients, 18.5% (n = 101/547) for painful-SFN patients, and 18.8% (n = 6/32) for painless-SFN patients). Of the variants detected, 70% were in SCN7A, SCN9A, SCN10A and SCN11A. The frequency of SCN9A and SCN11A variants was the highest in painful-SFN patients, SCN7A variants in painful-DPN patients, and SCN10A variants in painless-DPN patients. Our findings suggest that rare SCG genetic variants may contribute to the development of painful neuropathy. Genetic profiling and SCG variant identification should aid in a better understanding of the genetic variability in patients with painful and painless neuropathy, and may lead to better risk stratification and the development of more targeted and personalized pain treatments.

Muscle and Muscle-like Autoantigen Expression in Myasthenia Gravis Thymus: Possible Molecular Hint for Autosensitization.

The thymus is widely recognized as an immunological niche where autoimmunity against the acetylcholine receptor (AChR) develops in myasthenia gravis (MG) patients, who mostly present thymic hyperplasia and thymoma. Thymoma-associated MG is frequently characterized by autoantibodies to the muscular ryanodine receptor 1 (RYR1) and titin (TTN), along with anti-AChR antibodies. By real-time PCR, we analyzed muscle-CHRNA1, RYR1, and TTN-and muscle-like-NEFM, RYR3 and HSP60-autoantigen gene expression in MG thymuses with hyperplasia and thymoma, normal thymuses and non-MG thymomas, to check for molecular changes potentially leading to an altered antigen presentation and autoreactivity. We found that CHRNA1 (AChR-α subunit) and AIRE (autoimmune regulator) genes were expressed at lower levels in hyperplastic and thymoma MG compared to the control thymuses, and that the RYR1 and TTN levels were decreased in MG versus the non-MG thymomas. Genes encoding autoantigens that share epitopes with AChR-α (NEFM and HSP60), RYR1 (neuronal RYR3), and TTN (NEFM) were up-regulated in thymomas versus hyperplastic and control thymuses, with distinct molecular patterns across the thymoma histotypes that could be relevant for autoimmunity development. Our findings support the idea that altered muscle autoantigen expression, related with hyperplastic and neoplastic changes, may favor autosensitization in the MG thymus, and that molecular mimicry involving tumor-related muscle-like proteins may be a mechanism that makes thymoma prone to developing MG.

Integrative miRNA-mRNA profiling of human epidermis: unique signature of SCN9A painful neuropathy.

Personalized management of neuropathic pain is an unmet clinical need due to heterogeneity of the underlying aetiologies, incompletely understood pathophysiological mechanisms and limited efficacy of existing treatments. Recent studies on microRNA in pain preclinical models have begun to yield insights into pain-related mechanisms, identifying nociception-related species differences and pinpointing potential drug candidates. With the aim of bridging the translational gap towards the clinic, we generated a human pain-related integrative miRNA and mRNA molecular profile of the epidermis, the tissue hosting small nerve fibres, in a deeply phenotyped cohort of patients with sodium channel-related painful neuropathy not responding to currently available therapies. We identified four miRNAs strongly discriminating patients from healthy individuals, confirming their effect on differentially expressed gene targets driving peripheral sensory transduction, transmission, modulation and post-transcriptional modifications, with strong effects on gene targets including NEDD4. We identified a complex epidermal miRNA-mRNA network based on tissue-specific experimental data suggesting a cross-talk between epidermal cells and axons in neuropathy pain. Using immunofluorescence assay and confocal microscopy, we observed that Nav1.7 signal intensity in keratinocytes strongly inversely correlated with NEDD4 expression that was downregulated by miR-30 family, suggesting post-transcriptional fine tuning of pain-related protein expression. Our targeted molecular profiling advances the understanding of specific neuropathic pain fine signatures and may accelerate process towards personalized medicine in patients with neuropathic pain.

Large scale genome-wide association analyses identify novel genetic loci and mechanisms in hypertrophic cardiomyopathy.

Hypertrophic cardiomyopathy (HCM) is an important cause of morbidity and mortality with both monogenic and polygenic components. We here report results from the largest HCM genome-wide association study (GWAS) and multi-trait analysis (MTAG) including 5,900 HCM cases, 68,359 controls, and 36,083 UK Biobank (UKB) participants with cardiac magnetic resonance (CMR) imaging. We identified a total of 70 loci (50 novel) associated with HCM, and 62 loci (32 novel) associated with relevant left ventricular (LV) structural or functional traits. Amongst the common variant HCM loci, we identify a novel HCM disease gene, SVIL, which encodes the actin-binding protein supervillin, showing that rare truncating SVIL variants cause HCM. Mendelian randomization analyses support a causal role of increased LV contractility in both obstructive and non-obstructive forms of HCM, suggesting common disease mechanisms and anticipating shared response to therapy. Taken together, the findings significantly increase our understanding of the genetic basis and molecular mechanisms of HCM, with potential implications for disease management.

Semantic and right temporal variant of FTD: Next generation sequencing genetic analysis on a single-center cohort.

Semantic and right temporal variant of frontotemporal dementia (svFTD and rtvFTD) are rare clinical phenotypes in which, in most cases, the underlying pathology is TDP-43 proteinopathy. They are usually sporadic disorders, but recent evidences suggest a higher frequency of genetic mutations for the right temporal versus the semantic variant. However, the genetic basis of these forms is not clear. In this study we performed a genetic screening of a single-center cohort of svFTD and rtvFTD patients, aiming at identifying the associated genetic variants. A panel of 73 dementia candidate genes has been analyzed by NGS target sequencing including both causal and risk/modifier genes in 23 patients (15 svFTD and 8 rtvFTD) and 73 healthy age-matched controls. We first performed a single variant analysis considering rare variants and then a gene-based aggregation analysis to evaluate the cumulative effects of multiple rare variants in a single gene. We found 12 variants in nearly 40% of patients (9/23), described as pathogenic or classified as VUS/likely pathogenic. The overall rate was higher in svFTD than in rtvFTD. Three mutations were located in MAPT gene and single mutations in the following genes: SQSTM1, VCP, PSEN1, TBK1, OPTN, CHCHD10, PRKN, DCTN1. Our study revealed the presence of variants in genes involved in pathways relevant for the pathology, especially autophagy and inflammation. We suggest that molecular analysis should be performed in all svFTD and rtvFTD patients, to better understand the genotype-phenotype correlation and the pathogenetic mechanisms that could drive the clinical phenotypes in FTD.

The PINK1 p.Asn521Thr Variant Is Associated with Earlier Disease Onset in GRN/C9orf72 Frontotemporal Lobar Degeneration.

Genetic frontotemporal lobar degeneration (FTLD) is characterized by heterogeneous phenotypic expression, with a disease onset highly variable even in patients carrying the same mutation. Herein we investigated if variants in lysosomal genes modulate the age of onset both in FTLD due to GRN null mutations and C9orf72 expansion. In a total of 127 subjects (n = 74 GRN mutations and n = 53 C9orf72 expansion carriers), we performed targeted sequencing of the top 98 genes belonging to the lysosomal pathway, selected based on their high expression in multiple brain regions. We described an earlier disease onset in GRN/C9orf72 pedigrees in subjects carrying the p.Asn521Thr variant (rs1043424) in PTEN-induced kinase 1 (PINK1), a gene that is already known to be involved in neurodegenerative diseases. We found that: (i) the PINK1 rs1043424 C allele is significantly associated with the age of onset; (ii) every risk C allele increases hazard by 2.11%; (iii) the estimated median age of onset in homozygous risk allele carriers is 10-12 years earlier than heterozygous/wild type homozygous subjects. A replication study in GRN/C9orf72 negative FTLD patients confirmed that the rs1043424 C allele was associated with earlier disease onset (-5.5 years in CC versus A carriers). Understanding the potential mechanisms behind the observed modulating effect of the PINK1 gene in FTLD might prove critical for identifying biomarkers and/or designing drugs to modify the age of onset, especially in GRN/C9orf72-driven disease.

Peripheral Ion Channel Genes Screening in Painful Small Fiber Neuropathy.

Neuropathic pain is a characteristic feature of small fiber neuropathy (SFN), which in 18% of the cases is caused by genetic variants in voltage-gated sodium ion channels. In this study, we assessed the role of fifteen other ion channels in neuropathic pain. Patients with SFN (n = 414) were analyzed for ANO1, ANO3, HCN1, KCNA2, KCNA4, KCNK18, KCNN1, KCNQ3, KCNQ5, KCNS1, TRPA1, TRPM8, TRPV1, TRPV3 and TRPV4 variants by single-molecule molecular inversion probes-next-generation sequencing. These patients did not have genetic variants in SCN3A, SCN7A-SCN11A and SCN1B-SCN4B. In twenty patients (20/414, 4.8%), a potentially pathogenic heterozygous variant was identified in an ion-channel gene (ICG). Variants were present in seven genes, for two patients (0.5%) in ANO3, one (0.2%) in KCNK18, two (0.5%) in KCNQ3, seven (1.7%) in TRPA1, three (0.7%) in TRPM8, three (0.7%) in TRPV1 and two (0.5%) in TRPV3. Variants in the TRP genes were the most frequent (n = 15, 3.6%), partly in patients with high mean maximal pain scores VAS = 9.65 ± 0.7 (n = 4). Patients with ICG variants reported more severe pain compared to patients without such variants (VAS = 9.36 ± 0.72 vs. VAS = 7.47 ± 2.37). This cohort study identified ICG variants in neuropathic pain in SFN, complementing previous findings of ICG variants in diabetic neuropathy. These data show that ICG variants are central in neuropathic pain of different etiologies and provides promising gene candidates for future research.

Identification of risk loci for primary aldosteronism in genome-wide association studies.

Primary aldosteronism affects up to 10% of hypertensive patients and is responsible for treatment resistance and increased cardiovascular risk. Here we perform a genome-wide association study in a discovery cohort of 562 cases and 950 controls and identify three main loci on chromosomes 1, 13 and X; associations on chromosome 1 and 13 are replicated in a second cohort and confirmed by a meta-analysis involving 1162 cases and 3296 controls. The association on chromosome 13 is specific to men and stronger in bilateral adrenal hyperplasia than aldosterone producing adenoma. Candidate genes located within the two loci, CASZ1 and RXFP2, are expressed in human and mouse adrenals in different cell clusters. Their overexpression in adrenocortical cells suppresses mineralocorticoid output under basal and stimulated conditions, without affecting cortisol biosynthesis. Our study identifies the first risk loci for primary aldosteronism and highlights new mechanisms for the development of aldosterone excess.

Identification of a cytokine profile in serum and cerebrospinal fluid of pediatric and adult spinal muscular atrophy patients and its modulation upon nusinersen treatment.

Background and objectives: Multisystem involvement in spinal muscular atrophy (SMA) is gaining prominence since different therapeutic options are emerging, making the way for new SMA phenotypes and consequent challenges in clinical care. Defective immune organs have been found in preclinical models of SMA, suggesting an involvement of the immune system in the disease. However, the immune state in SMA patients has not been investigated so far. Here, we aimed to evaluate the innate and adaptive immunity pattern in SMA type 1 to type 3 patients, before and after nusinersen treatment. Methods: Twenty one pediatric SMA type 1, 2, and 3 patients and 12 adult SMA type 2 and 3 patients were included in this single-center retrospective study. A Bio-Plex Pro-Human Cytokine 13-plex Immunoassay was used to measure cytokines in serum and cerebrospinal fluid (CSF) of the study cohort before and after 6 months of therapy with nusinersen. Results: We detected a significant increase in IL-1ß, IL-4, IL-6, IL-10, IFN-γ, IL-17A, IL-22, IL-23, IL-31, and IL-33, in serum of pediatric and adult SMA patients at baseline, compared to pediatric reference ranges and to adult healthy controls. Pediatric patients showed also a significant increase in TNF-α and IL-17F levels at baseline. IL-4, IFN-γ, Il-22, IL-23, and IL-33 decreased in serum of pediatric SMA patients after 6 months of therapy when compared to baseline. A significant decrease in IL-4, IL-6, INF-γ, and IL-17A was detected in serum of adult SMA patients after treatment. CSF of both pediatric and adult SMA patients displayed detectable levels of all cytokines with no significant differences after 6 months of treatment with nusinersen. Notably, a higher baseline expression of IL-23 in serum correlated with a worse motor function outcome after treatment in pediatric patients. Moreover, after 6 months of treatment, patients presenting a higher IL-10 concentration in serum showed a better Hammersmith Functional Motor Scale Expanded (HFMSE) score. Discussion: Pediatric and adult SMA patients show an inflammatory signature in serum that is reduced upon SMN2 modulating treatment, and the presence of inflammatory mediators in CSF. Our findings enhance SMA knowledge with potential clinical and therapeutic implications.

Peripheral Ion Channel Gene Screening in Painful- and Painless-Diabetic Neuropathy.

Neuropathic pain is common in diabetic peripheral neuropathy (DN), probably caused by pathogenic ion channel gene variants. Therefore, we performed molecular inversion probes-next generation sequencing of 5 transient receptor potential cation channels, 8 potassium channels and 2 calcium-activated chloride channel genes in 222 painful- and 304 painless-DN patients. Twelve painful-DN (5.4%) patients showed potentially pathogenic variants (five nonsense/frameshift, seven missense, one out-of-frame deletion) in ANO3 (n = 3), HCN1 (n = 1), KCNK18 (n = 2), TRPA1 (n = 3), TRPM8 (n = 3) and TRPV4 (n = 1) and fourteen painless-DN patients (4.6%-three nonsense/frameshift, nine missense, one out-of-frame deletion) in ANO1 (n = 1), KCNK18 (n = 3), KCNQ3 (n = 1), TRPA1 (n = 2), TRPM8 (n = 1), TRPV1 (n = 3) and TRPV4 (n = 3). Missense variants were present in both conditions, presumably with loss- or gain-of-functions. KCNK18 nonsense/frameshift variants were found in painless/painful-DN, making a causal role in pain less likely. Surprisingly, premature stop-codons with likely nonsense-mediated RNA-decay were more frequent in painful-DN. Although limited in number, painful-DN patients with ion channel gene variants reported higher maximal pain during the night and day. Moreover, painful-DN patients with TRP variants had abnormal thermal thresholds and more severe pain during the night and day. Our results suggest a role of ion channel gene variants in neuropathic pain, but functional validation is required.

A multi-step genomic approach prioritized TBKBP1 gene as relevant for multiple sclerosis susceptibility.

BACKGROUND: Over 200 genetic loci have been associated with multiple sclerosis (MS) explaining ~ 50% of its heritability, suggesting that additional mechanisms may account for the "missing heritability" phenomenon. OBJECTIVE: To analyze a large cohort of Italian individuals to identify markers associated with MS with potential functional impact in the disease. METHODS: We studied 2571 MS and 3234 healthy controls (HC) of continental Italian origin. Discovery phase included a genome wide association study (1727 MS, 2258 HC), with SNPs selected according to their association in the Italian cohort only or in a meta-analysis of signals with a cohort of European ancestry (4088 MS, 7144 HC). Top associated loci were then tested in two Italian cohorts through array-based genotyping (903 MS, 884 HC) and pool-based target sequencing (588 MS, 408 HC). Finally, functional prioritization through conditional eQTL and mQTL has been performed. RESULTS: Top associated signals overlap with already known MS loci on chromosomes 3 and 17. Three SNPs (rs4267364, rs8070463, rs67919208), all involved in the regulation of TBKBP1, were prioritized to be functionally relevant. CONCLUSIONS: No evidence of novel signal of association with MS specific for the Italian continental population has been found; nevertheless, two MS loci seems to play a relevant role, raising the interest to further investigations for TBKBP1 gene.

Association of colorectal cancer with genetic and epigenetic variation in PEAR1-A population-based cohort study.

Platelet Endothelial Aggregation Receptor 1 (PEAR1) modulates angiogenesis and platelet contact-induced activation, which play a role in the pathogenesis of colorectal cancer. We therefore tested the association of incident colorectal cancer and genetic and epigenetic variability in PEAR1 among 2532 randomly recruited participants enrolled in the family-based Flemish Study on Environment, Genes and Health Outcomes (51.2% women; mean age 44.8 years). All underwent genotyping of rs12566888 located in intron 1 of the PEAR1 gene; in 926 participants, methylation at 16 CpG sites in the PEAR1 promoter was also assessed. Over 18.1 years (median), 49 colorectal cancers occurred, all in different pedigrees. While accounting for clustering of risk factors within families and adjusting for sex, age, body mass index, the total-to-HDL cholesterol ratio, serum creatinine, plasma glucose, smoking and drinking, use of antiplatelet and nonsteroidal anti-inflammatory drug, the hazard ratio of colorectal cancer contrasting minor-allele (T) carriers vs. major-allele (GG) homozygotes was 2.17 (95% confidence interval, 1.18-3.99; P = 0.013). Bootstrapped analyses, from which we randomly excluded from two to nine cancer cases, provided confirmatory results. In participants with methylation data, we applied partial least square discriminant analysis (PLS-DA) and identified two methylation sites associated with higher colorectal cancer risk and two with lower risk. In-silico analysis suggested that methylation of the PEAR1 promoter at these four sites might affect binding of transcription factors p53, PAX5, and E2F-1, thereby modulating gene expression. In conclusion, our findings suggest that genetic and epigenetic variation in PEAR1 modulates the risk of colorectal cancer in white Flemish. To what extent, environmental factors as exemplified by our methylation data, interact with genetic predisposition and modulate penetrance of colorectal cancer risk is unknown.

Dissecting the Polygenic Basis of Primary Hypertension: Identification of Key Pathway-Specific Components.

INTRODUCTION AND OBJECTIVES: Genome-wide association studies have identified a high number of genetic loci associated with hypertension suggesting the presence of an underlying polygenic architecture. In this study, we aimed to dissect the polygenic component of primary hypertension searching also for pathway-specific components. METHODS: The polygenic risk score (PRS) models, based on the UK biobank genetic signals for hypertension status, were obtained on a target Italian case/control cohort including 561 cases and 731 hyper-normal controls from HYPERGENES, and were then applied to an independent validation cohort composed by multi-countries European-based samples including 1,284 cases and 960 hyper-normal controls. RESULTS: The resulting genome-wide PRS was capable of stratifying the individuals for hypertension risk by comparing between individuals in the last PRS decile and the median decile: we observed an odds ratio (OR) of 3.62, CI = [2.01, 6.32] (P = 9.01E-07) and 3.22, 95% CI = [2.06, 5.10] (P = 6.47E-08) in the target and validation cohorts, respectively. The relatively high case/control ORs across PRS quantiles corroborates the presence of strong polygenic components which could be driven by an enrichment of risk alleles within the cases but also by potential enrichment of protective alleles in the old normotensive controls. Moreover, novel pathway-specific PRS revealed an enrichment of the polygenic signal attributable to specific biological pathways. Among those the most significantly associated with hypertension status was the calcium signaling pathway together with other mainly related such as the phosphatidylinositol/inositol phosphate pathways. CONCLUSIONS: The development of pathway-specific PRS could prioritize biological mechanisms, according to their contribution to the genetic susceptibility, whose regulations might be a potential pharmacological preventive target.

Congenital insensitivity to pain: a novel mutation affecting a U12-type intron causes multiple aberrant splicing of SCN9A.

ABSTRACT: Mutations in the alpha subunit of voltage-gated sodium channel 1.7 (NaV1.7), encoded by SCN9A gene, play an important role in the regulation of nociception and can lead to a wide range of clinical outcomes, ranging from extreme pain syndromes to congenital inability to experience pain. To expand the phenotypic and genotypic spectrum of SCN9A-related channelopathies, we describe the proband, a daughter born from consanguineous parents, who had pain insensitivity, diminished temperature sensation, foot burns, and severe loss of nociceptive nerve fibers in the epidermis. Next-generation sequencing of SCN9A (NM_002977.3) revealed a novel homozygous substitution (c.377+7T>G) in the donor splice site of intron 3. As the RNA functional testing is challenging, the in silico analysis is the first approach to predict possible alterations. In this case, the computational analysis was unable to identify the splicing consensus and could not provide any prediction for splicing defects. The affected intron indeed belongs to the U12 type, a family of introns characterised by noncanonical consensus at the splice sites, accounting only for 0.35% of all human introns, and is not included in most of the training sets for splicing prediction. A functional study on proband RNA showed different aberrant transcripts, where exon 3 was missing and an intron fragment was included. A quantification study using real-time polymerase chain reaction showed a significant reduction of the NaV1.7 canonical transcript. Collectively, these data widen the spectrum of SCN9A-related insensitivity to pain by describing a mutation causing NaV1.7 deficiency, underlying the nociceptor dysfunction, and highlight the importance of molecular investigation of U12 introns' mutations despite the silent prediction.

Correction: Evaluation of molecular inversion probe versus TruSeq® custom methods for targeted next-generation sequencing.

[This corrects the article DOI: 10.1371/journal.pone.0238467.].

Computational pipeline to probe NaV1.7 gain-of-function variants in neuropathic painful syndromes.

Applications of machine learning and graph theory techniques to neuroscience have witnessed an increased interest in the last decade due to the large data availability and unprecedented technology developments. Their employment to investigate the effect of mutational changes in genes encoding for proteins modulating the membrane of excitable cells, whose biological correlates are assessed at electrophysiological level, could provide useful predictive clues. We apply this concept to the analysis of variants in sodium channel NaV1.7 subunit found in patients with chronic painful syndromes, by the implementation of a dedicated computational pipeline empowering different and complementary techniques including homology modeling, network theory, and machine learning. By testing three templates of different origin and sequence identities, we provide an optimal condition for its use. Our findings reveal the usefulness of our computational pipeline in supporting the selection of candidates for cell electrophysiology assay and with potential clinical applications.